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Abstract 


The  dynamics  by  which  a  stationary  time  series  produces  its 
zero-crossings  and  higher  order  zero-crossings  sequentially  in 
time  is  studied,  illustrated,  and  applied  in  white  noise  testing. 
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1.  Introduction 

Consider  a  stationary  time  series  evolving  in  time  about 
level  zero.  As  the  time  series  threads  its  sample  path  about  this 
level  the  zero-crossing  rate  (the  number  of  zero-crossings  per 
unit  time)  converges  in  some  sense  as  the  series  length  increases. 
Similarly,  the  zero-crossing  rate  in  the  first  difference  of  the 
series  converges  too,  and  the  same  can  be  said  about  the  zero¬ 
crossing  rates  of  higher  order  differences  of  the  time  series.  In 
this  paper  we"  examine  the  dynamic  process  by  which  a  time  series 
and  its  higher  order  differences  produce  their  respective  zero¬ 
crossing  counts  sequentially  in  time  for  the  purpose  of  white 
noise  testing.  Because  a  zero-crossing  rate  as  a  function  of  the 
series  length  is  a  stochastic  process,  we  can  reiterate  by  saying 
that  we  shall  study  zero-crossing  rate  processes  obtained  from  the 
original  time  series  by  repeated  differencing.  Such  processes 
will  be  called  higher  order  crossings  rate  processes  or  simply 
HOC-processes .  The  present  work  is  an  extension  of  the  methods 
and  ideas  discussed  in  Kedem  (1987)  and  in  Kedem  and  Martin  (1987) 
where  HOC  plots  and  HOC-processes  are  introduced  as  graphical 
tools  in  time  series  model  identification.  Here  the  emphasis  is 
on  a  graphical  tool  that  is  more  powerful  than  HOC  plots  for  white 
noise  testing. 

There  is  a  difference  between  regular  higher  order  crossings 
(HOC)  and  HOC-processes.  HOC  are  zero-crossing  counts  in  filtered 
time  series  where  the  filtering  may  consist  of  any  combination  of 
linear  filters  applied  sequentially  or  repeatedly.  The  distinc¬ 
tion  lies  in  the  sample  size.  HOC  are  defined  for  fixed  record 


lengths  while  in  HOC-processes  the  record  lengths  increase.  In 
this  paper  by  HOC  we  always  mean  higher  order  crossings  obtained 
by  repeated  differencing.  General  properties  of  HOC  are  reviewed 
in  Kedem  ( 1986) . 

In  general  HOC  constitute  a  monotone  increasing  sequence. 

This  in  turns  implies  that  the  (k+l)'th  order  HOC-process  tends 
to  have  sample  paths  that  lie  above  the  sample  paths  of  the  k'th 
order  HOC-process,  a  fact  which  induces  a  certain  weak  convergence 
of  HOC-processes.  Thus  there  are  two  basic  modes  of  convergence 
associated  with  HOC-processes,  "horizontal"  that  depends  on  the 
series  size  and  "vertical"  that  depends  on  the  difference  order. 
Both  modes  of  convergence  summarize  the  oscillatory  behavior  of  a 
time  series.  The  initial  rate  of  the  vertical  convergence  is 
rather  fast  which  means  that  very  few  HOC-processes  are  needed  for 
our  purpose . 

The  first  part  of  the  paper  deals  with  basic  properties  of 
HOC  and  HOC-processes.  Section  4  discusses  in  some  detail  the 
difficulty  in  obtaining  the  exact  variance  of  HOC  and  suggests 
instead  a  useful  approximation.  Section  5  is  where  we  describe 
and  illustrate  our  test  procedure. 
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2.0  DEFINITION  OF  HOC  PROCESSES 


We  are  given  a  time  series  (Zt>,  t  =,*••,  -1 ,  0,  1 , 


with  covariance  function 


?( t,r)  =  Cov(Zt,Zt_r) ,  r , t  =  0,±  1,±  2, 


and  correlation  function 

p(t,r)  =  r (t,r) /r (t,0) ,  r , t  *  0,±  1,±  2, 


We  assume  the  process  is  stationary  with  zero  mean,  so  that 


=  E  Zt  >  E  Zt_J  =  E  Zt+j 


Mt,r)  =  E  {Zt  Zt_r)  -  rT 


rT  /  rQ. 


r  =  0,  ±1,  ±2, 
Gaussian. 


Finally  we  assume  that  the  process  {Zt>  is 


Let  {zj.0^}  s  { z^. }  be  the  0th-order  difference  process.  To 


form  the  higher  order  difference  processes,  we  begin  by  defining 


(1) 


the  first  order  difference  process  {Z^  )  such  that 


z(1)  =  z<0)  -  z(0)  *  V  z 

zt  zt  zt-l  zt' 


th 


In  general  for  the  k  -order  differences,  define 


z(k+l)  _  z|k)  _  z(k)  =  7(k+l)z^  k  a  0#1#  2 


(2.1) 


,(k) 


The  processes  (Z^.  }  are  referred  to  in  the  remainder  of  this 


paper  as  "Z-processes" . 


From  the  k-lth-order  difference  process  Z^k  ^  we  define 


the  kth-order  "clipped"  process  (x|k^  }  s’jch  that 


X^k)  =  I  (Zitk~1)  >  0),  t  -  0,  il,  .2, 


(2.2) 


s 


5 


where  !(•)  is  the  indicator  function.  The  binary-valued 
processes  {xjk^}  are  referred  to  as  "X-processes" .  We  further 
define 

d{.k)  =  I  [x*k)  *  x[kj]  .  (2.3 

(k.)  (  k-  1 ) 

Thus  d^  =  1  whenever  a  zero-crossing  occurs  between  zt_1 

( k-1 ) 

and  Z^.  ,  and  is  0  otherwise. 

The  total  number  of  zero-crossings  which  occur  between  t  =  1 
and  t  =  N  in  the  ( k-1 ) th-order  process  is  then 

VN  ’  4  +  d3  dN  '  |2'4 

The  processes  { N/N}  are  referred  to  as  HOC-processes  or 
simply  "D-processes" ,  with  k  being  the  order  and  N,  N  j  2,  the 
number  of  observations.  Thus  the  first  order  HOC-process  is 
*D1  '  N  =  2,3,* ••  .  More  precisely,  {D^  N/N}  18  a  sequence 

of  processes  with  parameter  N  indexed  by  k. 


3.0  PROPERTIES  OF  D-PROCESSES 


3.1  ACF  relationships 


Because  the  process  [Z 


(k+1) 


}  Is  obtained  from  simple 


differencing  of  the  process  {Z^.  '  },  there  is  a  direct  relation¬ 


ship  between  the  correlation  functions  of  the  Z-processes. 


/ =  covariance  function  for  {Z^k^ } 


=  correlation  function  for  {z|k^}. 


From  the  definition  of  {Z^.  }, 


z«k)  -Jo  (-D1  (*]  *w 


k  =  0,  l. 


so  that 


“jo  jo  ('1,1+3  (3  G9  E(Z‘-J  Z*-r-i  ' 


Rearranging, 


i- j=m 
Osisk 
Os  jsk 


We  know  that  for  positive  integers  a,  b,  n 


0  (n)  *  (3  (n0 . [n]  0  =  01 


and  we  observe  that 


■4.  0  0  -  0  0  ♦  [-5J . 0  00 


Os  is  k 
Os  Jsk 


(3-D 


,  -  ,  -  J.  -  s  ‘  .  /,  «\  •'  -  .  -  .  ■*.  ’  .**« 

>  v\-  -  .-.v.  -  v.-- 
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■  ty  (9  -  U-S-x]  (ty 

-  K-J 


(o)  (ky 


Therefore , 


"  ■  P) '™  P.) 


i-.i- 


and  we  can  obtain  an  expression  for  p^  ;  setting  r  *  1 


(k)  _ 


'2k 

k-1 

+pl( 

'2k1  f2k  ’ 
k  k-2 

» _ i 

)-P2( 

'2k  ' 
k-1 

+ 

'2k 

k- 

(f) 

-2Pl 

f2k  I 

.k-lJ 

+  2P2 

p.) 

“2P3 

'2k  'J 
k-3  J 

)+•••+(-!)  p 


+• • -+(-1)^2 p. 


(3-2) 


For  r  2:  1 ,  we  can  obtain  from  the  basic  relationships 

„ ( k+1 )  _ ( k+1 )  . 

(k+1)  *  '  zt  t-r  ' 


so  that 


E  {  (z{k+1))2  } 


E  {  (Z<k)  -  Z<k>)  (Z<k>  -  2^^) 
E  {  (z|k)-  zJkJ)2  } 


2  v<k>  -  v<k)  -  y(k> 

rr  Ar  +  1  r t-1 

,  v(k)  (k) 

*  *0  *  *1 


2  p<k+1>  (1  -  p<k>)  =  2p(k)  -  P(k)  -  P(ki 

"  '  1  '  ^ r  r+l  r-1 


P«J>  =  -  P'k|  +  2  Pi">  -  2p<k+1>  II  -  p'k)>  (3.3) 

r+l  r-l  r  r  1 


3.2  Expected  Values  for  D-Processes 


When  the  underlying  proces  (Zt)  is  Gaussian,  the  expected 
number  of  zero  crossings,  E  N,  can  be  evaluated  directly  as  a 
function  of  p*k  1 ^ ,  the  first  correlation  coefficient  of 
We  observe  that 


and  that 


,(k)  _  ,x(k)_y(k).2 
t  “  (Xt  Xt-lJ 


E  d{k)  =  E  { (x|k)-xjkJ)2} 


=  E  {(x[k))2  +  (X[kj)2  -  2  x|k)  x[k|  ) 


( k)  .  2 


■  (k)  v  ^  „/v(k) 
t-1 

since  the  x-processes  are  0/1 -valued. 


=  E  (X^.  '  }  +  E{X|_')  -  2  E  (xj"'  X^J  > 


(k)  v(k) 
‘t 


(3.4) 


The  assumption  that  {Zt)  is  zero-mean  Gaussian  implies  that 
(k) 

(Zt  >  is  also  zero-mean  Gaussian,  k  ■  1,  2,***,  since  from 
(2.1)  the  higher  difference  processes  are  linear  combinations  of 
the  0-order  process .  Then 

E  (X[k))  =  P  (x|k)  =  1} 

=  P  (z£k_1)  *  0) 

=*  1/2 

and 


E  <X<k)  x[k>>  =  P  {  ( z|k  1]  i  0)  and  (Z<*~1}  s  0)) 


=  1/4  +  (l/27r)  sin  1  pjk  1] 


(3.5) 


Thus  we  have 


The  relationship  (3.7)  together  with  (3.2)  and  (3.3)  indica¬ 
tes  that  the  expected  zer  -j-crossing  rates  are  uniquely  determined 
by  the  correlation  structure  of  the  underlying  Z^0^-  process. 


3.3  Convergence  of  HOC-processes 

It  is  instructive  to  consider  the  graphs  of  some  HOC- 
processes.  Figures  3. 1-3. 3  show  the  graphs  of  several  HOC- 
processes  obtained  from  both  stalonary  as  well  as  nonstationary 
time  '■■'ries.  The  first  striking  observation  is  that  the  HOC- 
processes  seem  to  converge  monotonically  to  a  limiting  HOC-process 
and  that  the  convergence  is  rather  fast.  This  is  the  vertical 
convergence  mentioned  earlier,  and  one  that  has  been  observed  in 
numerous  cases  as  already  noted  in  Kedem  and  Martin  (1987).  The 
horizontal  convergence  is  observed  in  the  figures  at  least  in  the 
stationary  cases.  We  shall  now  discuss  these  observations  more 
formally  below. 

A  basic  fact  which  underlines  the  vertical  convergence  is 


that 


(3.8) 
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From  this  inequality  we  can  show  that  for  strictly  stationary 
processes  the  expected  HOC  constitute  a  monotone  sequence. 


UN  HU'  PROCESS*  DI-D8  ) 


Figure  3.1.  HOC-processes  obtained  from  two  white  noise  time 
series . 


a 

Theorem  3.3.1  obviously  covers  stationary  Gaussian  processes,  but 
it  is  only  a  first  step  in  the  explanation  of  the  vertical  con¬ 
vergence.  For  a  general  result  we  may  appeal  to  the  theory  of 
point  processes  as  outlined  in  Kallenberg  (1976). 

It  is  convenient  to  restrict  attention  to  the  set  [0,®) 
equipped  with  the  Borel  field  ©.  Let  IB  be  the  class  of  subsets 
consisting  o~  all  bounded  sets  in  ®.  The  set  of  all  locally 
finite  measures  that  are  nonnegative  and  integer  valued  is  denoted 
by  N .  Assume  that  the  process  {Z^}  is  defined  over  the 


probability  space  (Q,A,P)  and  introduce  in  M  the  a-algebra  7 


generated  by  the  mappings  //  — *  ti(  B),  ft  e  M ,  Be®.  By  a  point 

process  on  [0,®)  we  mean  any  measurable  mapping  from  (Q,j4,P)  into 

With  an  obvious  extension  of  the  previous  notation,  let 

D^CB),  B  e  (B  ,  be  the  number  of  zero-crossings  generated  by 
k—  1 

V  in  B.  Then  for  each  realizations  of  (zt)  »  Dk(*) 

defines  a  measure  in  M  and  so,  for  each  fixed  k,  is  a 

point  process.  The  vertical  convergence  can  be  explained  with  the 
help  >1  the  sequence  of  point  processes  {D^(*)>.  We  need  the 
fol  1  jwing  combinatorial  lemma. 

Lemma  3.3.1.  For  any  process  (Z^)  there  exists  a  tQ,  l5tQ«N 
and  tQ  depends  on  j,  j=l,2,--*,  such  that  for  every  N 


Z  d{J>  ,  £ 

_ 4  L  A.  _  4 


(3.9) 


t*l 
t*  t„ 


with  probability  one. 


Proof  :  Fix  j.  Suppose  d| • +d^^  =  0;  then  we  are  done. 

So  suppose  d j ^ +• • • +d^-^  >  0.  Then  there  exists  at  least  one 
symbol  change  in  (x|^  )  ,  l<t5N.  Note  that  in  general  we  have  the 
implication 

(X^J  *  X*J))  c  +  =  (xj^). 

It  follows  that  there  exists  t 1  such  that 

x[J+1>  =  x<J>  . 

1  *1 


Therefore  each  symbol  change  in  (X|J/)  produces  at  least  one 


symbol  change  in  (X 


( J  +  l) 


}  for  tj<tsN.  On  the  other  hand  by 


going  from  j  to  j+1  we  may  lose,  by  (3.8),  at  most  one  zero¬ 
crossing  for  lstst^  It  follows  that  there  exists  a  tQ  which 
depends  on  j  such  that  (3.9)  holds.  a 


The  lemma  says  that  in  the  binary  sequence  (x|^+1^)  ,  lstsN, 
there  are  N-l  locations  that  give  at  least  as  many  symbol  changes 
as  the  corresponding  locations  in  (X^^ }  ,  IstsN,  while  no 
information  is  available  about  the  remaining  location  which  we 
denote  by  tQ.  With  the  help  of  the  lemma  we  can  now  prove  the 
weak  convergence  of  the  sequence  of  point  processes  (D.  (•)}. 

Theorem  3.3.2.  Let  {Zt}  b*  a  strictly  stationary  process  and 
let  D(*)  be  the  point  process  corresponding  to  the  limiting  HOC- 


a 


Di  s  Di+i  surelv- 

/  lr  ) 

By  strict  stationarity  and  since  d^  is  binary 

,(,™ , ,) ,  , ,) 

Therefore 


or 


(Di +  k’ 1 »] 5  p(D>Ui -  d^+1)  8  0 

pk„  -  °]  *  pK«,k  -  °] 


and  so  the  sequence  P(Dk  N  “  °)  converges  monotonically  as  k  — ►  ® 
for  any  N.  By  appealing  to  Theorem  4.7  in  Kallenberg  (1976)  it 
follows  that  D^( • )  converges  in  distribution  to  D(*)*  □ 

For  the  horizontal  convergence  we  require  ergodlcity. 


Theorem  3.3.3.  Let  (Zt)  be  an  ergodic  zero-mean  stationary 
Gaussian  process.  Then  for  every  k. 


11m 

N-*oo 


D 


k,  N 


N 


1  -1 (k-1 ) , 

-  cos  (pj  ) 


a .  s . 


Proof :  Apply  (3.7). 


□ 


This  explains  the  convergence  to  straight  horizontal  lines 


observed  in  the  figures.  It  follows  that  for  Gaussian  white  noise 


the  straight  lines  are  at  levels  -  cos  1(-k/(k+l)),  k=0,l,2, 

n 


When  both  the  order  and  the  size  increase  we  have  convergence 
towards  the  highest  possible  frequency  in  the  spectrum  if  the  later 


exists.  More  precisely  we  have.  Define 

D. 


d3(N)  s  -2^  . 
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Also,  the  normalized  spectral  measure  of  is  given  by 


yj(do>)  = 


j^sln  -  dF(h>) 

Sn  [sin  |  xj  dF ( X 


where  F  Is  the  spectral  distribution  function  of  { Z ^ }  -  Let  &u 
be  the  point  mass  at  u. 

Theorem  3.3.4.  Let  {Z^.}®  be  a  2ero  raean  stationary  Gaussian 

* 

process,  and  suppose  o  s  n  is  the  highest  frequency  in  the 

_  * 

spectral  support.  If  lim  D,(N)  s  w  /n  with  probability  one, 

N-*»  3 

j  =  1 , 2  ,  • • •  ,  then 

* 

lim  lim  D.(N)  *  —  a.s. 
j-too  N-»®  3  n 

Proof.  From  Kedem  and  Slud  (  1982)  *  I  6  *  +  *  6  *,  j  -♦  ®. 

J  Z  —(a)  Z  (t) 

But  the  Gaussian  assumption  implies  that 

(n 

cos(wEDj(N))  *  J  cos (u ) v j ( dw )  — ►  cos(w  ),  j  — ►  ®. 


Therefore  ED ^ ( N )  — >  <o  /n ,  j  — ►  ®,  uniformly  in  N. 

By  strict  stationarity  we  can  define  the  a.s.  limit 

c  .  s  lim  D, (N) . 

3  N-tao  3 

* 

Then  by  Lemma  1  c .  <  c.  a.s.,  and  so  lim  c  .  ®  c  <  w  /n  a.s. 

3  J  +  1  j-#x>  3 

Thus  by  bounded  convergence  E(Cj)  — *  E(c),  j  — ►  ®.  Strict 
stationarity  and  Fatou's  Lemma  yield 

E(c.)  2  lim  sup  ED.(N). 

3  N-too  3 

_  * 

Since  ECj  — ►  Ec  and  EDj(N)  — ►  w  /n , 
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E(  c  )  2  cj  /n  . 


*  * 

But  c  s  q  /n .  Therefore  c  =  <*>  /n  a.s.  □ 

* 

Corollary  3.3.1.  If  u  =  n  then 

11m  lim  D  (N)  =  1. 
j-*oo  N-wo  ^ 

We  close  this  section  by  noting  that  Theorem  3.3.2  which 
addresses  the  vertical  convergence  is  somewhat  more  general  than 
the  Higher  Order  Crossings  Thcjrem  proved  in  Kedem  and  Slud  (1982) 
and  which  says  that  for  a  strictly  stationary  process  with 
spectral  distribution  function  F  supported  at  n,  as  k  — ►  oo, 

rx(k)1  •••  010101  •••  w.p.  1/2 

t  J  •  •  •  101010  •  •  •  w.p.  1/2 

Theorem  3.3.2  does  not  presuppose  that  {Zt)  possesses  moments  at 
all,  and  strictly  speaking,  it  is  not  a  spectral  result  but  can  be 
viewed  as  such  in  a  generalized  sense.  It  is  important  to  observe 
that  a  process  need  not  possess  moments  of  any  order  and  yet  its 
HOC-processes  possess  moments  of  all  orders. 
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4.0  ESTIMATION  OF  VARIANCE  OF  D-PROCESSES 


4.1  Direct  Computation 

In  order  to  use  zero-crossing  statistics  for  testing,  one 

needs  at  least  variance  estimates  for  the  D.  These  estimates 

k ,  N 

are  in  general  functions  of  fourth  order  orthant  probabilities  in 
which  are  not  easily  obtainable. 

In  addition  calculation  of  the  fourth  order  orthants  leads  to 
a  difficult  computational  problem  since  the  error  appears  to 
increase  rapidly  with  N. 

In  what  follows  we  discuss  in  some  detail  the  problem  of 
approximating  the  variance  of  ^ .  Fortunately,  it  is  possible 

to  obtain  a  very  reasonable  approximation  under  a  certain 


condition. 


4.2  Simulation 


As  an  alternative  to  direct  computation  of  the  variance, 
approximations  to  the  variance  of  D,  „  can  be  obtained  from 
computer  simulations.  The  computer  is  programmed  to  generate 
realizations  of  a  specified  process  from  which  counts  of  the 
higher  order  zero-crossings  are  made.  Counts  are  obtained  for  a 
number  of  realizations,  and  the  sample  mean  and  variance  of 
these  statistics  are  calculated. 

In  our  study,  we  have  used  computer  simulation  results  as  a 
comparison  for  variance  estimates  obtained  by  analytical  methods. 
We  have  implemented  these  simulations  both  on  a  large  mainframe 
computer  (UNI VAC  1108)  using  FORTRAN,  and  on  an  IBM  personal 
computer  using  TurboPASCAL.  We  note  that  while  the  execution  time 
is  considerably  faster  on  the  mainframe,  running  time  on  a  PC  is 
not  prohibitively  slow  (approximately  20  minutes  for  100  itera¬ 
tions  with  1000  observations  and  six  levels  of  higher-order 
crossings)  and  offers  some  advantages  of  flexibility  and  lower 
cost . 

Some  examples  of  simulation  results  are  shown  in  Tables 
4- 1 A , B , C  and  4-2A,B,C.  The  results  agree  with  an  earlier 
simulation  reported  in  Kedem  (1987) 


TABLE  4-1 A 


Simulation  Results  for  D 

J.  /  w 

White  Noise  Process 
(100  iterations) 


N 

®  °1,N 

HDi.n))1/2 

10 

4.71 

1 . 54 

20 

9.96 

2.03 

50 

24.93 

3.46 

100 

49 .61 

4.92 

200 

99.46 

7.84 

500 

249.76 

10.16 

1000 

499.15 

15.28 

TABLE  4-1B 

Simulation  Results  for  D 

AR(1)  Process  $*0.5 
(100  iterations) 


[varK«))1/2 


10 

3.65 

1 . 44 

20 

6 . 88 

2.55 

50 

16.52 

3.53 

100 

32.77 

4 . 89 

200 

67.43 

6.88 

500 

167 . 56 

11.64 

000 

336.89 

15.91 

TABLE  4- 1C 


Simulation  Results  for  N 

AR( 1 )  Process  0=0.2 
(100  Iterations) 


N 

E  D1,N 

[Var[°l,N]J 

10 

4.40 

1.63 

20 

8.58 

2 . 10 

50 

22 . 56 

3.60 

100 

43.51 

5.08 

200 

86.78 

7.05 

500 

217.84 

10.11 

1000 

437.83 

15.76 

jjVr, 


TABLE  4-2A 

Simulation  Results  for 

White  Noise  Process 
(N=1000,  100  iterations) 


500.16 

16.23 

668 . 00 

12.05 

732 . 26 

11.86 

768.49 

11.24 

793 . 20 

10.93 

811.74 

10.64 

TABLE  4-2B 

Simulation  Results  for 

AR(l)  Process  0=0.5 
(N=1000,  100  iterations) 


335.16 

16.58 

578.87 

13.54 

683 . 52 

13.94 

739. 10 

12.00 

773.83 

11.64 

797.95 

10.64 

4.3  Approximation  Assuming  Negligible  Comulants 

For  any  stationary  (Z^.)  (not  necessarily  Gaussian)  with 
known  mean  ,  define 


Ch  =  C-h  =  T^h  E  (Zt-/0  (Z1 


h  =  0 ,  1  ,  •  •  •  ,  T-  1 


which  is  an  unbiased  estimate  of  /  . 

n 

It  is  shown  in  Anderson  (1971,  p.452)  that 


T-h-1 


(T-h)  Var  C  =  £  (l  -  III  *  I 


r=- ( T-h- 1 ) 


T-h-J  L/  r  '  r  +  h  '  r-h 


+  #(h. -r,h-r) ] 


where  h  is  the  cumulant  function 


*  (  h ,  r  ,  s  )  =  E  {(Zt  -  //)  (Zt+h-/0  (Zt  +  r-/i)  (Zt+s-/i)) 

'  t/h  /  r_s  +  /  r  /'h_s  +  t  s  7-h_r] 
t  =  • • • ,  1,  0,  1, 

When  (Z^)  is  Gaussian  x  vanishes. 

We  observe  that 


D)'N  '  1=2  ^ 


~  V  f  V  (  1  )  Y  (  1  )  1 

'  A  ri  “  Xi-ll 


2  N 


-&K’) 


N 

i  x|U  -  1/2]  |xj  1/2  +  (N-1J/2 


(N-l)  2C* 


(  1  ) 


N-l 

2 


X  '  (  1  ) 

where  is  the  function  (4.1)  for  the  Xv  -process.  Thus, 

D  1  M  1  f"  y  (  1  ) 

Var  [^rrj  =  4  Var  [C1  J  •  {4-3) 

If  we  can  assume  that  the  cumulants  are  negligible  for  the 
processes  we  are  interested  in,  then  (4.2)  is  easy  to  evaluate  and 
can  be  used  to  give  estimates  of  the  variance  of  N  by  using 

(4.3)  . 

We  calculated  estimates  of  the  variance  of  D,  for 

1  ,  N 

Gaussian  White  noise  and  for  Gaussian  AR(1)  processes  with  para¬ 
meters  0.5  and  0.2.  The  results  are  shown  in  Table  4-3.  For 
white  Gaussian  noise  the  cumulants  are  0,  so  (4.2)  works  exactly 
for  this  process.  However,  the  estimates  for  the  AR(1)  processes 
appear  to  overestimate  the  variance  compared  to  the  simulation 
results  in  Table  4-2B.C,  and  the  differences  are  larger  as  the 
AR(1)  process  moves  away  from  white  noise.  Thus  we  conclude  that 
this  method  does  not  provide  good  variance  estimates  for  processes 
other  than  white  noise.  In  other  words,  the  cumulant  function 
contribution  to  (4.2)  is  appreciable. 


4.4  Approximation  for  White  Noise 


Kedem  and  Reed  (1986)  showed  that  for  an  m-dependent 
stationary  Gaussain  process  (Z^)  with  «  in  the  spectral  support 
and  for  fixed  N 


1  im 

k-*v 


var(DkV 


(N-l )A  ( 1 - A  |k)  ) 


1 


(4.4) 


where 


A jk)  =  P(xJk)  =  1 


=  1/2  +  1 /»  sin  1  pjk 


k=l , 2 , * • ■ 


Thus  approximation  to  the  variance  of  ^  is  provided  by 

Var(Dk  )  =  (N-l)  a  <k)  (1-A<k))  (4.5) 

This  method  is  incorporated  in  the  method  using  first  order 
Markov  approximation  (Section  4.5).  We  may  view  (4.5)  as  0 1  th 
order  approximation. 


4.5  Approximation  assuming  first  order  Markov  process. 

The  difficulties  in  the  computation  of  Var(Dk  )  can  to  a 

( k ) 

large  extent  be  passed  by  assuming  that  the  binary  process  (d^  )  , 

t = 0  ,  *■  1  ,  •  •  •  ,  is  a  Markov  chain  for  each  fixed  k.  This  idea  has 
been  suggested  and  implemented  in  Kedem  (1987)  where  it  is  shown 
that  the  resulting  approximation  is  surprisingly  good.  The  method 
is  summarized  in  the  following  algorithm. 

Given  ’  ’  ,f\‘  Pjc+1  >  we  can  compute  E  ^  exactly  and 

approximate  Var  ( ^)  by  following  steps  1-7. 


Step  1  .  Obtain 


from  (3.2) . 


Step  2.  Obtain  from  (3.3). 


Step  3.  Compute 


.  k  1  1  .  -1  (k- 1  •  ,  „ 

A.  =  -  +  -  sin  />.  '  ,  j  =  1 , 2 

J  2  71  j  J 


Step  4 .  Compute 


Step  5 .  Compute 


l-*Jk)  ,  q(k) 


= 


l-2A<k)  +  i‘k) 
2( 1-A  Jk)  ) 


Step  6.  The  exact  expected  value  of  N  is  given  by 

E  Dk  N  =  (N-l)p(k» 

Step  7.  The  variance  of  Dk  N  is  approximated  by 


Var(Dk  N)  =  (N-l)p(k)q(k)  + 


2p<kVk)(,(k)-p<k>) 


|<N-1»“Vk.N} 


where 


vk.N  '  ^^TkT, 

The  variance  expression  (4.6)  can  be  recognized  as  the 


variance  of 


,(k)  ,(k) 

d2  +---+  dN 


(4.6 


under  the  assumption  that  (d|  )  is  a  two  state  stationary 


(  k  i  (  k  ) 

Markov  chain  with  parameters  p  and  r 


Observe  that  (4.6! 


is  (4.5)  plus  a  "correction"  term.  Tables  4.4A  -  4 . 4D  show  some 
computations  using  formula  (4.6).  These  results  agree  closely 
with  simulation  results.  Subroutine  DVAR  attached  at  the  end, 


uses  the  above  algorithm  in  the  computation  of  (4.6). 


TABLE  4-4 A 


Results  of  First  Order  Markov 
Approximation  for  White  Noise 
(N  =  1000) 


k 


E  D 


k ,  N 


1  /2 


1 

499 . 50 

15.80 

2 

666 . 00 

13.14 

3 

731.55 

12.17 

4 

769 . 18 

11.58 

5 

794 .37 

11.16 

6 

812.76 

10.82 

TABLE  4 - 4  B 


Results  of  First  Order  Markov 
Approximation  for  AR  (  1  )  </>  =  0.5 
(N  =  1000) 


3  1 


k 

E  D.  M 

Var 

°k,  N 

J 

1/2 

1 

333.00 

15.7b 

2 

579 . 85 

13.99 

3 

684 . 69 

12.74 

4 

729 . 67 

12.00 

5 

773.82 

11.49 

6 

797 . 44 

11.10 

TABLE  4-4C 


Results  of  First  Order 
Approximation  for  AR(1) 
(N  =  1000) 

Markov 

0=0.2 

k 

E  Vn 

IMvn 

1 

435 . 47 

15.82 

2 

630 . 36 

13.41 

3 

709 . 85 

12.39 

4 

754 . 45 

11.77 

5 

783 . 61 

11.33 

6 

804 .48 

10 . 97 

5.  Testing  for  white  noise:  an  application  of  HOC  processes 

Finally,  after  discussing  the  dynamic  behavior  of  HOC- 
processes  in  regard  to  order  and  time,  we  are  ready  to  apply  this 
interesting  dynamics  in  tests  appropriate  for  time  series  models. 

We  concentrate  on  testing  for  white  noise  but  the  technique  can  be 
readily  extended  to  more  general  tests  regarding  autoregressive- 
moving  average  processes. 

As  noted  earlier,  HOC-processes  summarize  the  oscillatory 
characteristics  of  a  process.  Moreover,  it  is  apparent  from 
Figures,  3. 1-3. 2  that  for  stationary  time  series,  HOC-processes 
produce  sample  paths  which  settle  about  straight  lines  even  for 
moderate  sample  sizes,  and  that  the  initial  rate  of  the  vertical 
convergence  is  rather  fast  which  is  a  general  property  of  higher 
order  crossings.  Thus,  in  this  sense,  it  is  possible  to  characte¬ 
rize  a  process  by  a  few  HOC-processes  observed  between,  say, 

N  =  500  to  N  =  1000.  Experience  shows,  and  this  will  be  illustrated 
again  below,  that  as  few  as  six  HOC-processes  need  to  be  considered. 
Based  on  these  observations,  the  proposed  test  for  white  noise  is 
to  observe  whether  the  first  six  HOC-processes  produce  sample  paths 
which  indeed  fall  within  prescribed  probability  bounds  derived 
under  the  hypothesis  of  white  noise.  The  hypothesis  of  white  noise 
is  rejected  as  soon  as  at  least  one  of  the  six  paths  cross  its 
probability  bounds.  Such  a  scheme  is  akin  to  control  processes 
and  can  be  used  in  real  time  monitoring  of  man  made  systems. 

Observe  that  under  the  Markov  assumption,  N  has  an 
asymptotic  normal  distribution  so  that  it  lies  in 


LiLV.WlC'^V.V.I^lVlV 


with  95%  approximately,  for  each  fixed  k  and  large  N.  However 


as  N  varies  we  must  obtain  a  simultaneous  bound  for  the  HOC  path 


{ j^/N )  ,  -  N  •  N2 .  In  this  work  we  take  N  ^  =  500,  N2=1000, 


l<k<6.  The  bounds  are  determined  from 


c/Var  (  ^ )  ,  500-  N-  1000 


for  each  k,  l'k'6,  where  (Var(D^  ^ ) } 


is  given  by  (4.6) 


assuming  white  noise.  The  constant  c  controls  the  overall  power 


and  the  a-level  of  the  test.  When  c=1.96  the  test  is  very 


powerful  against  many  alternatives  but  also  has  a  high  a-level. 


In  order  to  lower  the  level  we  widen  the  bounds  by  choosing  larger 


values  for  c,  a  procedure  that  leads  to  more  conservative  bounds. 


5.1  Power  considerations  and  examples 


Simulations  were  first  run  to  estimate  the  power  of  the  test 


for  white  noise  as  a  function  of  c,  the  test  width  coefficient  in 


(5,2).  Three  processes  were  simulated:  white  noise  and  two  AR ( 1 


processes  with  parameters  0.05  and  0.1.  The  resulting  HOC  process 


paths  were  monitored  using  the  white  noise  limits  obtained  from 


(4.6)  and  (5.2).  Figures  5.  1-5.3  show  the  percentage  of  cases 


failing  the  test  for  the  three  simulations. 


From  Figure  5.1  it  appears  that  a  c- value  close  to  3  corres- 


pons  to  a  level  around  5%.  In  fact,  we  obtained  an  a-level  of 


0.06  for  C  =  2 . 8 .  The  power  corresponding  to  c  =  2.8  is  about  20% 


and  40%  for  0=0.05  and  0=0.1,  respectively.  This  is  shown  in 


Figures  5.2,  5.3.  The  curves  shown  in  the  three  figures  were 


obtained  by  a  least -squares  fit 


,%  *.  *.  v  *. 


Table  5.1  gives  the  results  of  some  power  calculations 
obtained  from  100  simulations  for  each  alternative  from  ARMA(1,1). 
The  value  of  c  is  2.8.  The  further  the  alternative  is  from 
white  noise  the  higher  is  the  power  as  expected. 

Examples  which  illustrate  the  use  of  the  test  procedure  are 
shown  in  Figures  5. 4-5. 6.  For  non-white  noise  in  many  cases 
higher  order  HOC-processes  resemble  those  from  white  noise  and  the 
distinction  is  made  by  low  order  HOC-processes .  Occasionally  it 
is  the  D6-path  that  captures  the  difference  while  the  lower  order 
cases  do  not . 
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Table  5.1.  Power  for  the  white  noise  test  against  alternatives 
from  Z^.  =  0Zt  t-  u^_  -  Ou  where  u^  are  independent  N(0,1) 
random  variables. 


/.  ■*,  V,  »*, 


Figure  5.4  Hoc-processes  from  white  noise.  The  hypothesis  of 
white  noise  is  accepted.  Each  HOC-path  is  within  the  white  noise 
1 imi ts . 
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Figure  5 . 6  Hoc-processes  from  ARMA(1,1),  <fi=t>  =  0.1,  and  white  noise 
limits.  The  hypothesis  of  white  noise  is  rejected. 
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5.2  Identification  of  general  models. 

There  is  no  reason  why  the  above  procedure  cannot  be  extended 
to  test  for  non-white  models,  because  once  a  model  is  completely 
specified  and  its  correlation  structure  is  known  as  a  result, 

(4.6)  and  (5.2)  can  be  evaluated  3nd  the  limits  can  be  obtained  as 
before.  Our  experimental  results  indicate  that  to  maintain  a 
reasonable  test  level  the  value  of  2.8  for  c  in  (5.2)  is  still  a 
good  choice.  Table  5.2  gives  some  power  calculation  for  testing 
the  hypothesis  that  a  process  is  AR(1)  with  parameter  <P  ~  0.1, 
where  the  alternatives  are  generated  by  ARMA ( 1 , 1 )  with  parameters 
0  ,0  as  above  . 
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0 . 51 

0 . 30 

0 . 00 

0 . 97 

0 . 05 

0 . 50 

1 . 00 

Table  5.2.  Power  of  the  test  for  testing  the  hypothesis  that  the 
process  is  AR { 1 )  with  0  =0.1.  The  alternatives  are  generated  by 
different  ARMA(1,1)  processes. 


We  have  outlined  a  test  procedure  for  time  series  models 
based  on  HOC-processes.  The  test  can  be  used  in  online  situations, 
for  example  when  monitoring  the  oscillation  of  signatures  obtained 
from  man  made  systems  such  as  engines.  As  was  indicated,  if 
higher  a-levels  can  be  tolerated,  the  power  against  different 
alternatives  can  increase  appreciably.  This  can  be  achieved  by 
lowering  c  and  thus  rendering  the  limits  less  conservative. 

The  ideas  presented  in  this  paper  can  be  extended  to  cover 
nonstationary  time  series  as  well  since  the  vertical  convergence 
is  quite  a  general  phenomenon  and  seems  to  hold  across  a  wide 
range  of  processes . 

The  present  paper  is  a  generalization,  an  application  of  HOC 
to  one  based  on  HOC-processes.  For  each  fixed  N  the  sequence 
{D^  ,  k  =  1 , 2 , • • • ,  is  the  sequence  of  higher  order  crossings. 

It  is  known  that  the  rate  of  increase  in  HOC,  and  hence  also  in 
HOC-processes,  is  a  useful  feature  in  discriminating  between 
processes  as  is  exemplified  in  Figures  3. 1-3.3.  Tests  based  on 
{ D i  n,*,,,D6  for  fixed  N  are  discussed  in  Kedem  (1987).  It 
seems  that  no  real  advantage  is  gained  by  increasing  the  order  k 
in  the  present  use  of  HOC-processes  amidst  stationary  time  series. 
On  the  other  hand  by  increasing  N  we  can  obtain  a  more  powerful 
test.  This  pertains  particularly  to  borderline  cases  that  need 
"more  time"  to  get  out  of  the  probability  bounds  as  they  become 
narrower  and  narrower. 

Finally,  the  method  described  herein  impinges  on  sequential 
analysis.  This  connection  will  be  pursued  elsewhere. 


C  SUBROUTINE  DVAR  COMPUTES  E ( D( K ) ) , VAR ( D ( K ) ) , K= 1 , . . . , 0 
C  GIVEN  THE  FIRST  7  AUTOCORRELATIONS  OF  A  ZERO  MEAN 
C  STATIONARY  GAUSSIAN  TIME  SERIES. 

C  R: INPUT. VECTOR  CONTAINING  FIRST  7  AUTOCORRELATIONS. 

C  N : INPUT . NET  LENGTH . LENGTH  OF  TIME  SERIES  MINUS  5  (DIFFERENCES) 

C  RDsOUTPUT. VECTOR  CONTAINING  R D ( K ) ,  K  =  0 , ,6. 

C  RD ( K )  IS  1  *  ST  AUTOCOR.  IN  THE  K ' TH  DIFFERENCE . 

C  RD(K)  IS  COMPUTED  FROM  R ( 1 ),..., R ( K+ 1 ) 

C  ED:OUTPUT. VECTOR  CONTAINING  EXPECTED  HOC  K=  1  , . .  . , 6 

C  VD: OUTPUT. VECTOR  CONTAINING  VARIANCES  OF  HOC  K=I,...,6. 

★ 

SUBROUTINE  DVAR ( R,N, RD, ED, VD ) 

REAL  R ( 7 ) , E D ( 6 ) , VD ( 0 ) , RD (0:6) 

R  D  (  0  )  =  R  (  1  ) 

RD(  1  )=(-l.+2.*R(  1  )  -  R  ( .1 )  )/(  2 .  *  (  I .  -  R  (  1  )  )  ) 

RD(  2)  =  (-  4.+7.  *R(  1  )-4.  *R(2)+R(  3)  )/(6.-8.  *R(  1  )+2.  *K(  2)  ) 

RD(  3 )  =  (  - 1 5 . +  26 . *R (  1 )  - 1 6 . *R ( 2 ) +6 .  *R ( 3 ) -R  (  4 )  )/ 

1  (  20.-30. *R( 1 )+12.*R(2)-2.*R(3)  ) 

RD ( 4 )  =  ( -56 . +98 . *R ( 1 ) -64 . *R ( 2 ) +29 . *R ( 3 ) -8 . *R ( 4 ) + R ( 5  )  )/ 

1  (70.-1 12. *R( 1)+50.*R(2)-16.*R(  3)+2.*R(4)  ) 

RD(  5)  =  (-210.+372.  *R(  1  ) -255 .  *R  (  2  )  + 1  30 .  *R  (  3  )  -4(> .  *K  (  4  )  +  1 0 .  *H  (  5  ) 

1  -R(b) )/( 252.-420. *R( 1 ) + 240 . *R ( 2 ) -90 . * R ( 3 ) +20 . *R ( 4 ) -2 . *R ( 5 ) ) 

RD ( 6 )  =  ( -  792 . + 1 4 1 9 . *R ( 1 ) - 1 01 2 . *R ( 2 ) +561 . *  R (  3 ) -2  32 . *R ( 4 ) +67 . *R ( 5 
1  -  1 2 . *R ( 6 ) +R ( 7 ) )/( 924.-1584. *R( 1 ) +990 . *K ( 2 ) -440 . *R ( 3 )  + 

1  132.*R(4)-24.*R(5)+2.*R(6) ) 

DO  1  1=0,5 
R 1 =RD( I ) 

R2=-l.+2. *RD (I)— 2. *RD( 1  +  1 )*( 1. -RD ( I  )  ) 

ALAM 1=0. 5+ASIN ( R 1 )/ACOS(-l . ) 

ALAM2=0 . 5+ AS IN ( R2 ) / ACOS ( - 1 . ) 

P=1 . -ALAM 1 
Q=ALAM 1 

PLAM= (1.-2. * ALAM 1 + ALA M2 ) / ( 2 . * ( 1 . -ALAM 1 )  ) 

V 1 = ( X- 1 ) *  P  *  o 

V2  =  2  .  *P*0* ( DLAM-P )/(  1 . -DLAM ) 

V3=Q* ( l.-(  ( DLAM-P ) /O ) *  * ( N- 1  )  )/(  1 .-DLAM) 

V4=V2* (N-1-V3 ) 

VD ( I + 1 ) =V1+V4 
ED ( 1+1 )  =  (N-1  )  *P 

i  conti  r:uF. 

RETURN 

END 
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